Astronomy & Astrophysics manuscript no. 9309manu 


© ESO 2008 


February 2, 2008 





Letter to the Editor 

Detection of diffuse interstellar bands in M31 

M. A. Cordiner 1 , N. L. J. Cox 2 , C. Trundle 1 , C. J. Evans 3 , 1. Hunter 1 , N. Przybilla 4 , F. Bresolin 5 , and F. Salama 6 

1 Astrophysics Research Centre, School of Mathematics and Physics, Queen's University, Belfast, BT7 INN, U.K. 

2 Herschel Science Centre, European Space Astronomy Centre, ESA, P.O.Box 78, E-28691 Villanueva de la Canada, Madrid, Spain. 

3 UK ATC, Royal Observatory Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, U.K. 

4 Dr. Remeis-Sternwarte Bamberg, Sternwartstr. 7, D-96049 Bamberg, Germany 

Institute for Astronomy of the University of Hawaii, 2680 Woodlawn Drive, 96822 Honolulu, HI, U.S.A. 
6 Space Science Division, NASA Ames Research Center, Mail Stop 245-6, Moffett Field, California 94035, U.S.A. 

Received December 21st 2007 / Accepted January 23rd 2008 

ABSTRACT 

Aims. We investigate the diffuse interstellar band (DIB) spectrum in the interstellar medium of M31. 

Methods. The DEIMOS spectrograph of the W. M. Keck observatory was used to make optical spectroscopic observations of two 
supergiant stars, MAG 63885 and MAG 708 17, in the vicinity of the OB78 association in M3 1 where the metallicity is approximately 
equal to solar. 

Results. The /U5780, 5797, 6203, 6283 and 6613 DIBs are detected in both sightlines at velocities matching the M31 interstellar Nai 
absorption. The spectra are classified and interstellar reddenings are derived for both stars. Diffuse interstellar band (DIB) equivalent 
widths and radial velocities are presented. 

Conclusions. The spectrum of DIBs observed in M31 towards MAG 63885 is found to be similar to that observed in the Milky 
Way. Towards MAG 708 17 the DIB equivalent widths per unit reddening are about three times the Galactic average. Compared to 
observations elsewhere in the Universe, relative to reddening the M31 ISM in the vicinity of OB78 is apparently a highly favourable 
environment for the formation of DIB carriers. 
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extinction 



1. Introduction 

Today, more than 300 diffuse interstellar bands (DIBs) are 
known but the carriers have remained unidentified since their 
discovery almost 100 years ago. It is debated whether the DIB 
carriers arise from the dust, the gas, or the large-m olecule com - 
ponent of the interstellar medium (see the review bv lSarrel2006l) . 
The s ubstructure present in many of the DIB pro files indi- 
cates (IS aire et al. 1 119951: lEhrenfreund & Foind fl996t) that they 
are caused by large gas-phase molecules. The interaction be- 
tween gas-phase species (atoms and molecules), UV radiation 
and dust grains may be crucial in the formation of large (> 50 
atom) molecules in space which are likely candidates for the car- 
riers o f the DIBs (see e.g. ISalama et al.lll996l: iRuiterkamp et al.l 
120051) . 

In interstellar environments distinct from those found in the 
Milky Way (MW), previous research on the relationships be- 
tween atoms, molecules, dust and DIBs has focused on the Large 
and Small Magellanic Clouds (e.g. Ehrenfreund et al.l 120021 : 
ICox et al.ll2rj06Ti2007t IWeltv et al.ll2006l) . The behaviour of the 
DIB carriers was analysed with respect to the higher gas-to- 
dust ratios, lower metallicities, lower Ry and stronger interstellar 
radiation fields of these environments. Beyond the Magellanic 
Clouds, studies are sparse, co nfined to sightlines probed by suffi- 
ciently bright supernovae (e. g.lRichll987l:ISqllerman et al.l2005h 



or background quasars (e.g. lYork et al.ll2006l) 



M3 1 presents the opportunity to study the effects on the DIB 
carriers of the unique chemical and physical conditions found 
in this Local Group galaxy. This letter presents spectra of two 
supergiant stars in M3 1 and we report the first unambiguous de- 
tection of DIBs in M3 1 . The observed DIB properties are dis- 
cussed in relation to those in other galaxies and the interstellar 
conditions of M3 1 . 



2. Observations 

Seve nty-two bright stars i n and around the M31 OB 78 associa- 
tion (Ivan den Berghfl964l) were observed in November 2003 us- 
ing the Keck DEIMOS spectrograph. Two angles of the 1200G 
grating were used to cover the region from approximately 3500 
to 9000 A with a resolving power R - 3300. The total expo- 
sure time was 2.25 hours in the blue region and 1.5 hours in 
the red, during which the seeing was 0.5 - 0.8". 1 Wavelength 
calibration accuracy was confirmed using the Na D sky emis- 
sion lines whose central wavelengths were found to be accurate 
to within +5 km s _1 . Reduced spectra were Doppler-corrected 
to the LSR frame. Table [T] shows the co-ordinates and photom- 
etry of MAG 63885 and MAG 70817. These targets were se- 
lected for analysis in this letter due to their prominent DIBs. The 



Send offprint requests to: Martin Cordiner, e-mail: 
m . cordiner@qub .ac.uk 



1 Initially the red and blue spectra were both reduced using the 
spec2d pipeline. However, the wavelength scale of the blue spectra 
(~ 3500 - 6000 A) was found to be erroneous so these were re-reduced 
manually in IRAF using the twodspec routines. 
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Table 1. Observed M31 stars. Target numbers from lMagnier et all ( |1992|) . Co-ordinates and photometry from lMassev et al.ld2006l) . 
Deri ved stellar spect ral classifications and radial velocities (v») are shown. M31 Eg-v values were calculated using intrinsic colours 
from lJohnsonl ( 1 19661) and have been corrected for Galactic foreground reddening (see Sect. 13.4b . 
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Fig. 1. Observed continuum-normalised Nai spectra (his- 
tograms) with best-fitting models overlaid (red curves). Velocity 
scale at the top is relative to the Na D[ rest wavelength. The M3 1 
and Galactic (MW) absorption components are labelled. 



rest of the DEIMOS data will be presented in a future article 
(ICordiner et al.ll2008l) . 

3. Analysis and results 

3.1. Interstellar sodium D 

The observed Na D lines (see Figure [U show two clearly sepa- 
rated absorption components due to interstellar gas at velocities 
consistent with the Galactic and M31 ISM (around radial ve- 
locities of and -55 km s" 1 respectivel y). The spectra were 
modelled using vapid dHowarth et a l. 2002) and each absorption 
component was accurately reproduced with a Gaussian interstel- 
lar cloud model. The fits are plotted in Figure [TJ and the mean 
radial velocities and equivalent widths (EWs) of the M31 Di ab- 
sorption components are given in Table [2] 

3.2. Diffuse interstellar bands 

We searched the spectra for all DIBs with central depths greater 
than .05 in Figure 6 2 of the DIB survey bv lJenniskens & Desertl 
(I1994I) . Sc aled to the reddening of our targets, this central depth 
limit corresponds to ~ Icr of the spectral Poisson noise. The 
/U5780, 5797, 6203, 6283 and 6613 DIBs were detected at M31 
velocities in the spectra of MAG 63885 and MAG 70817. -16196 
was tentatively detected towards MAG70817. The /U4501, 
4726, 5705, 5849, 6269, 6376, 6379, 6445, 6532, 6660, 6993 
and 8026 DIBs were too weak to be detected. Telluric absorp- 
tion line contamination prevented analysis of the /L16886, 6919, 
7224 and 7334 DIBs. 

To simultaneously measure the DIB radial velocities and 
equivalent widths, Galactic DIB templates were fitted to the ob- 
served spectra using a non-linear least-squares algorithm. The 
Galactic templates wer e derived from h igh resolution, high S/N 
spectra of /3 l Sco (see ICordinerll2006l) . These were shifted to 



the interstellar rest frame and convolved with a Gaussian in- 
terstellar cloud model then with the instrumental spectral PSF 
and rebinned to the wavelength scale of the DEIMOS spectra. 
This technique minimises statistical errors in the measurement of 
the equivalent widths and velocities of weak DIBs but assumes 
that the intrinsic M3 1 DIB profiles and rest wavelengths closely 
match those of /?' Sco. If this is not the case, unknown system- 
atic errors could occur, but we find no evidence to suggest that 
M3 1 DIB profiles differ at all from those of /?' Sco. Observed 
DIB spectra and fitted profiles are shown in Figure 12 

For the /W5780, 6283 and 6196 DIBs the radial velocities 
and equivalent widths were allowed to vary in the fits. A telluric 
absorption component (derived from a hig h-resolution spec trum 
of the unreddened fast-rotator a Gru (see ICordiner||2006l) . de- 
graded to R = 3300), was also included for /16283 to account 
for the prominent band-head of atmospheric O2 between 6278 
and 6286 A. 3 For the AA5797, 6196 and 6203 DIBs only the 
equivalent widths were allowed to vary while the DIB radial ve- 
locities were fixed at the M31 Nai radial velocities. In all cases, 
there is excellent agreement between the fitted and observed DIB 
profiles. DIB radial velocities and EWs are shown in Table [2] 
The radial velocities of all measured DIBs closely match the 
Nai radial velocities, consistent with previous high-resolution 
studies of the Galacti c an d extragalactic ISM (s uch as those by 
ISollerman et ap2005l and ICox et al.ll2006ll2007l) . 

Statistical equivalent width and DIB velocity error estimates 
were derived by replicating the spectra 100 times with randomly 
generated Poisson noise added to each. The DIB model param- 
eters were refitted for each replication and individual parameter 
errors were obtained from the +68th percentiles of the resulting 
parameter ranges. 

We also found evidence for a broad absorption feature to- 
wards MAG 63885 centred near 4430 A, identified as a blend of 
M31 and Milky Way ,14428 DIBs. This feature was modelled 
as a sum of two Lorentzian components with FWHM 17.25 A 
(see lSnowll2002l) at v = and v = -550 km s" 1 . The best-fitting 
model yields 1057 mA for the M31 component and 236 mA for 
the MW foreground component. However, this result is uncertain 
due to stellar line contamination and the relatively poor signal- 
to-noise in this spectral region. 

3.3. Stellar spectral types 

Spectral types and stellar radial velocities for the target stars are 
presented in Table [TJ Analysis of luminous B-type supergiants 
in M31 finds chemical a bundances comparab le to those in the 
solar neighbourhood (e.g. iTrundle et al.ll2002T) . so Galactic stan- 
dards are appropriate. The spect ra were classified with ref erence 
to the Galactic standards from lEvans & Howarfhl (120031) . with 
lumino sity classes assign ed on the basis of the Hy equivalent 
widths (lEvans et al.ll2004l) . 



2 A synthetic average of their observed DIB spectra (normalised to 
E B - V = 1.28). 



3 Our DEIMOS observations contained no telluric standards by 
which to perform a conventional telluric division. 



Cordiner, Cox, Trundle, Evans, Hunter, Przybilla, Bresolin, & Salama: Detection of diffuse interstellar bands in M31 



3 



^5780 LSR velocity (km s ) 
-1500 -1000 -500 500 1000 1500 




5750 5760 5770 5780 5790 5800 5810 
LSR wavelength (A) 

X6203 LSR velocity (km s" 1 ) 
-1500 -1250 -1000 -750 -500 -250 250 500 750 



|5Sco 



MW6196 MW6203 




6170 6175 6180 6185 6190 6195 6200 6205 6210 6215 6220 
LSR wavelength (A) 




6260 6265 6270 6275 6280 6285 6290 6295 
LSR wavelength (A) 



)u6613 LSR velocity (km s ) 
-1500 -1250 -1000 -750 -500 -250 250 500 750 



3 °- 95 - 63885 „ n 

| 0.85 ^JU^^f 



6580 6585 6590 6595 6600 6605 6610 6615 6620 6625 6630 6635 
LSR wavelength (A) 

Fig. 2. Continuum-normalised spectra of the observed DIBs 
(black histograms). Galactic (MW) and M31 DIBs are labelled. 
Prominent stellar lines are marked with asterisks. Fitted model 
DIB profiles are overlaid (red curves). Velocity scales at the top 
are relative to the Galacti c DIB rest wavelengths published by 
iGalazutdinov et all (l2000h . /3 l Sco Galactic reference DIB tem- 
plates (plotted with arbitrary intensity scaling factors), are shown 
at the top of each figure, with additional telluric C»2-band spec- 
trum (dashed line) for the /I6283 region. 



Table 2. Summary of best-fitting DIB equivalent widths 
(EW /mA) and LSR radial velocities (v/kmf 1 ) for MAG 63885 
and MAG 708 17. / denotes values held fixed during fitting. 
Statistical (68th percentile) DIB velocity errors are less than 
about +10 kms -1 . The M31 Nai radial velocities and D; line 
EWs are also shown, including an estimated contribution due 
to stellar photospheric Na i in square brackets. EWs for /3 l Sco 
(Eb-v = 0.22) and Galactic mean EW/Eb-v da ta (normalised 
to E B -y - 0-14) are also given (da ta from iHer bi g l I1993L 
iThorburn et ai]|2003llMegier et alJl2005l and lCordiiieril2006l) . 
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Stellar radial velocities are the mean results from man- 
ual measurements of the cores of the Ha Balmer line and the 
Paschen lines. Statistical errors on the measured radial veloci- 
ties are less than +5 km s . 

We searched Galactic stellar spectra of the same spectral 
types as our targets for the presence of lines overlapping the de- 
tected DIBs. No significant contamination of the AA6196, 6203, 
6283 or 6613 DIBs is expected. The AA5780 and 5797 may suf- 
fer contamination of up to about 5 and 2 mA respectively as a 
result of overlapping lines of Fe i. 

3.4. Foreground gas and dust 

LAB H i data (iKalberla et al]l2005l) for the nearest survey point 
in the direction of our sightlines shows N(Hi) = 5.91 x 10 20 
cmT 2 over the velocity ra nge of Galactic gas (from -300 to 
100 km s~'). Equation 7 of iBurstein & Heilesl dl978l) is used to 
calculate the foreground reddening for lines of sight at latitudes 
away from the Galactic plane which yields Eb-v = 0.06 mag to- 
wards our t argets, identical to the value given by the foreground 
dust map of lSchlegel etail d 19981) . 

The Galactic foreground Na i column densities measured to- 
wards MAG 63885 and MAG 708 17 are 3.7 x 10 12 and 1.9 x 10 12 
cm -2 resp ectively, whic h correspond to reddenings of 0.12 and 
0.08 mag dHobbsll 19741) . However, there may be contamination 
of the Galactic Na i profiles near v = due to sky-line subtrac- 
tion residuals. Using a foreground reddening of 0.06, corrected 
M31 Eb-v values are given in TableQ] 

4. Discussion 

DIB equivale nt widths in the Milky Way correlate with Eb-v 
dHerbidl995l) . Given the signal-to-noise and the limited redden- 
ing of our target spectra, only the strongest known Galactic DIBs 
were detected in M31. Other M31 DIBs (including those listed 
in Section [3~!2l are expected in our spectra around or below the 
limit of detectability. 

The DIB equivalent widths observed in M31 toward 
MAG 63885 are consistent with those typically observed in 
Galactic sightlines with the same reddening, as highlighted by 
comparison with the similarly-reddened /J 1 Sco (see Table |2J. 
Towards MAG 708 17 the observed DIBs are stronger per unit 
Eb-v by a factor of about 2-5 compared with Galactic-mean 
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Fig. 3. /15780 equivalent widths vs. redd ening for differen t 
galaxie s Data for the LMC and SMC (IWelty et al.1 l2006h . 
DLAs (lYork et al.1 l2006t lEllison et al.l |2008j) and NGC 1448 
dSollerman et al.l 20051) are s hown. The dashed line shows a lin- 
ear fit to Galactic d ata (fromlHerbisll993L iThorburn et"aT]|2003L 
iMegieret al.ll2005l and lCordinerll2006l) . 



data. Figure|3]shows that the ^5780 EW/E B -v for MAG 63885 is 
among the larger values observed elsewhere in the Universe and 
for MAG 708 17 /15780 EW/E B _ V is significantly greater than 
any other sightline plotted. 

Detailed information is sparse, but M31 appears be similar to 
the Milky Way in terms of the average interste llar gas-to-dust 
ratio (Ne dialkov et al.ll2000l iBresolin et al.ll2002l) and the metal- 
licity in the vicinity of our targets ([Trundle et al.ll2002l) . There 
is, however, evidence for a difference in the properties of the in- 
terstellar dust grains and the interstellar UV radiation field. A 
different star f ormation history and low rate of star formation 
( 1/ 1 0th of MW; IWalterbos & Braunl 19941) has been measured in 
M3 1 , and the surface r adiation flux was fo und to be poor in UV 
dCesarskvet al.l 1 19981; iPagani et al.lll999l) . The fact that strong 
DIBs are observed under these conditions may be in contradic- 
tion to the hypothesis that UV radiation is required for the p ro- 
duction of the carriers (see lHerblil l 1995HKendall et al.ll2002l for 
example). However, MAG 63885 and MAG 70817 are near to 
the OB78 association where the abundance of early-type stars 
may result in a strong interstellar UV radiation field. 

The M31 2175 A UV extinction bump has been found to 
be weak and narrow dBianchi et alJll996tlHutchings et al.|[l 992) 
and a peculiar extinction law was observed by iMassev et al.l 
d 1995b in the anomalously low average colour-excess ratio 
Eu-bIEb-v- A ccepting curre nt theories of dust grain extinction 
(see review by lDraing|2003l) . the evidence is consistent with a 
different distribution of dust grain sizes and compositions com- 
pared to the Milky Way average. In partic ular, there may be a 
lack of small graphitic dust grains (see also lXu&Heioulll994l) . 
In that case, the observation of strong DIBs implies that the car- 
riers are not closely associated with the small grains believed to 
be responsible for shape of the UV extinction curve. 

5. Conclusion 

The ^5780, 5797, 6283, 6203 and 6613 DIBs were detected to- 
wards MAG 63885 and MAG 708 17 at velocities corresponding 
(within 25 km s~') to the mean M31 interstellar Nai absorption. 
/16196 and /14430 were also tentatively detected. 

The M31 DIB spectrum towards MAG 63885 is consistent 
with that observed in the Galaxy. Towards MAG 708 17 the DIBs 
are about a factor of three stronger per unit reddening than the 
Galactic average. The high DIB strengths might be related to 



differences in the M3 1 interstellar UV extinction curve and ra- 
diation field compared to the Galaxy. Further studies will be re- 
quired to determine if the two sightlines examined here are rep- 
resentative of the general M3 1 DIB behaviour. 

A planned instrument for the Hubble Space Telescope, the 
'Cosmic Origins Spectrograph' will be able to provide essential 
information on the UV properties of M3 1 sightlines to further 
examine these possibilities. However, the connection between 
DIBs in M31 and the atomic and molecular content {e.g. Na, 
K, Ca, CH + , CH, CN and C 2 ) of the diffuse medium can only 
be properly addressed with high resolution, high signal-to-noise 
optical spectroscopic observations. 
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